Universal Breakdown of Elasticity at the Onset of Material Failure 
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We show that, in the athermal quasi-static deformation of amorphous materials, the onset of 
failure is accompanied by universal scalings associated with a divergence of elastic constants. A 
normal mode analysis of the non-affine elastic displacement field allows us to clarify its relation to 
the zero-frequency mode at the onset of failure and to the crack-like pattern which results from the 
subsequent relaxation of energy. 



Experiments on nanoindentation of metallic glasses , 
on granular materials and on foams Q, demonstrate 
that at very low temperature and strain rates, the mi- 
crostructural mechanisms of deformation involve highly 
intermittent stress fluctuations. These fluctuations can 
be accessed in molecular dynamics simulations, but are 
best characterized numerically via "exact" implementa- 
tion of a-thermal quasi-static deformation: alternating 
elementary steps of affine deformation with energy relax- 
ation permits one to constrain the system to reside 
in a local energy minimum (inherent structure) at all 
times. As illustrated in figure ^ macroscopic stress fluc- 
tuations arise from a series of reversible (elastic) branches 
corresponding to deformation-induced changes of local 
minima. These branches are interrupted by sudden ir- 
reversible (plastic) events which occur when the inher- 
ent structure annihilates during a collision with a saddle 
point. H These transitions constitute the most elemen- 
tary mechanism of deformation and failure for disordered 
materials at low temperature. 

Using this quasi-static protocol, recent studies of both 




0.29497 0.29498 0.29499 
7 



FIG. 1: Stress (top) and shear modulus (bottom) for a small 
strain interval about a strain of .3. Left: fixed strain steps of 
size 10" 4 ; Right: convergence to the yield point (circled on 
the left) with decreasing strain step. 



elasticity and plasticity could identify important 
properties of elasto-plastic behavior which arise solely 
from the geometrical structure of the potential energy 
landscape. Tanguy et al Q have observed that, follow- 
ing reversible (elastic) changes of the inherent structures, 
molecules undergo large scale non-affine displacements. 
They have shown these non-affine displacements to be 
related to the breakdown of classical elasticity at small 
scales and to quantitative differences between measured 
Lame constants and their Born approximation. Malan- 
dro and Lacks have shown that the destabilization of 
a minimum occurs through shear-induced collision with 
a saddle. At the collision, a single normal mode sees 
its eigenvalue going to zero. Building on this work, we 
studied the irreversible (plastic) event following the dis- 
appearance of an inherent structure: subsequent mate- 
rial deformation in search of a new minimum involves 
non-local displacement fields -in the likeness of nascent 
cracks- controlled by long-range elastic interactions. 

Several molecular displacement fields thus appear to 
be closely related to the geometrical structure of the 
potential energy landscape: (i) non-affine displacements 
along elastic branches, (ii) the single normal mode con- 
trolling the annihilation of an inherent structure, and (iii) 
the overall deformation occurring during an irreversible 
event. In order to piece together a complete picture of 
elasto-plasticity at the nanoscale, we need to understand 
the relation between these different fields and ask how 
elastic behavior breaks down at the onset of failure. It is 
thus a study of incipient plasticity -that is the onset of 
irreversible deformation- that we wish to perform. Here, 
the structural disorder is expected to control the onset 
of failure: this situation is somehow opposite to homo- 
geneous defect nucleation in crystals, [S| where failure is 
controlled by Hill's continuum condition. 0] 

We base our approach on exact microscopic expressions 
for the non-affine corrections to elasticity in disordered 
solids, [TnL Hll | which have entirely been overlooked in 
recent works. Here, we put such analytical tractations 
in perspective with the recent numerical developments. 
We derive an exact formulation for the non-affine dis- 
placement fields, and construct a normal mode decom- 
position therefrom. This analytical framework permits 
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us to evidence that the lowest frequency normal mode 
dominates the non-affine elastic displacement field close 
to a plastic transition. We then show that at any plastic 
transition point, energy stress, and the vibration fre- 
quencies display a singular, universal, behavior associ- 
ated with a divergence of the elastic constants. The nor- 
mal mode analysis of the subsequent cascade shows that 
the low frequency modes active at incipient plasticity are 
superseded by long-range elastic interactions in the latter 
stages of an irreversible event. 

We consider in this work a molecular system in a peri- 
odic cell. The geometry of the cell is determined by the 
matrix h whose columns are the Bravais vectors. ^2, El 
The affine deformation of the cell between configurations 
ho and h is characterized by the Green-St Venant strain 
tensor, e = | ((Tig 1 ) T .h T .h.h^ 1 — l), which governs the 

elongation of a vector x — * x, as x 2 = x 2 + 2x T .e.x. 
As the energy functional generally depends only on the 
set of interparticle distances, it can be parameterized as 
U({?i\, e) where {fi} are the positions of the particles in 

a reference cell. ^3 Q] Varying, e for fixed {r^} corre- 
sponds to an affine displacement of the molecules in real 
space. 

To start, let's contemplate more closely the athermal, 
quasi-static algorithm. Deformation, e(j), is enforced by 
moving the Bravais axes of the periodic cell; 7 is intro- 
duced as rescaled coordinate which measures the amount 
of deformation from some reference state. In practice 
5(7) corresponds to either pure shear or pure compres- 
sion. Formally, the limit ho = h (or 7 — > 0) is often ap- 
propriate to define stresses and elastic constant around 
a (possibly stressed) reference state. Once a choice of 
ho is made, the algorithm tracks in the reference cell 
a trajectory {^}(7)i which is implicitly defined by de- 
manding that the system remain in mechanical equilib- 
rium: [io,E3i 
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Starting at mechanical equilibrium at 7 = 0, {^(7)} is a 
continuous function of 7 on some interval [0,7 C ]. At j c , 
the local minimum collides with a saddle point. [f| 

An equation of motion for f = {^(7)} is obtained 
by derivation of |QJ with respect to 7. Denoting, H — 
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This relation holds for any 7 £ [0,7c]- In the limit 
h — > hg, H. is the Dynamical Matrix. To invert 7i, trans- 
lation modes must be eliminated by fixing the position of 
a molecule, df / d7 is a rescaled "velocity" of molecules in 
quasi-static deformation. It defines the direction (in tan- 
gent space) of the non-affine displacement field observed 





FIG. 2: Left: The force response to simple shear, H, at a 
strain configuration, 7 = 0.2945, or 7 C — 7 ~ 10~ 4 . Right: 

dr 



The non-affine velocity (or "displacement") field, 4^ for the 



same state as shown on the left. 



by Tanguy et al and can be directly evaluated by solv- 
ing equation J5J without resorting to quadruple precision 
minimization, [6j 

Here, we illustrate these ideas with numerical simula- 
tions of a two-dimensional bidisperse mixture of parti- 
cles interacting through a shifted Lennard- Jones poten- 
tial. Particle sizes r$ — r^sin f^/sin ^ and a number 

ratio Nl/Ns = 1+ 4 V ^ are used to prevent crystallization; 
the simulation cell is 50 xtl in length. We have also 
performed simulations on Hertzian spheres to check that 
it yielded results consistent with those presented here. 
Typical patterns of the fields H and dr/d7 in (steady) 
simple shear deformation are shown figure the appar- 
ent small scale randomness of the vector S is in sharp 
contrast with the large vortex-like structures displayed 
by the non-affine "velocity" field dr/d7- To understand 
the randomness of H, note that Ej = dFi/dj is the force 
response on molecule i after an elementary affine defor- 
mation of the system: it only depends on the configu- 
ration of the molecules with which molecule i interacts, 
hence is an e-dependent measure of the local disorder of 
molecular configurations. We checked that spatial cor- 
relations decay very fast in the field H: in the following 
discussion, the short-range randomness of the field H al- 
lows us to interpret it as noise. 

An analytical expression for the bulk elastic constants 
derives along similar lines. [lollll| The first derivative of 
the potential with respect to the components of e defines 
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the thermodynamic stress, t: t a p — de - — de 
The total derivative indicates derivation while preserving 
mechanical equilibrium, the second equality results from 
equation (JTJ , and Vo is the volume of the simulation cell. 
The second (total) derivative of the energy gives the elas- 
tic constants, 1 1 111 
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We recognize the first term as being the Born approxima- 
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tion C^p"^ . The second term accounts for the non-affinc 

corrections, and reads: C a p xcr = — ^a a p.H.~ 1 .a xa . 
Similarly, the second derivatives of the energy, following 
any generic deformation £(7), can be written as: 



d 2 U d 2 U 



d-y 2 dj 2 
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For an isotropic material, the elastic constants can be 
written: C a( 3 X(7 = XS a pS xa + ^,(S ax Si3a + S acr S/3 X ), which 
define the Lame constants, A and /x. In order to estimate 
these constants, it is not necessary to evaluate all the 
components of the tensor S = (S a/ g), but only two of 

its projections H, e.g. for pure shear and pure compres- 
sion and use equation (0J . In equation (0} the correction 
to the Born term is negative definite: quantities such as 
the shear modulus, /x, or the compressibility, K = A + /x, 
are necessarily smaller than the Born term, while this is 
not necessarily true of A = K — /x alone as it does not, 
by itself, correspond to any realizable mode of deforma- 
tion. This is consistent with the numerical observations 
by Tanguy et al 6] in Lennard- Jones systems. 

Next, we perform a normal mode analysis of the fields 
H. Denoting the eigenvectors of the Dynamical Ma- 
trix (normal modes), and A p the associated eigenvalues, 
the vector H can be decomposed as: H = £p ^p; with 

£ p = H.\& p . (If H is a random field, the variables £ p 
are random.) From this decomposition, expressions can 
be obtained for the non-afnne direction and for the non- 
affinc contribution to elasticity: 
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We now concentrate on the behavior of the shear mod- 
ulus at incipient plasticity, as shown in figure ^ an d [3] 
Malandro and Lacks have shown numerically that at 
the onset of a plastic event a single eigenfrequency goes 
to zero. HI We denote **(7) the first non-zero normal 
mode; in two dimensions, it is the third in the spectrum. 
Close to failure (7 — > -f c ), A* (7) — > 0, hence **(7 C ) must 
dominate the non-afnne direction dr/d7. This is true if 
and only if the quantity = V&*(7).E(7) does not 

vanish at the yield point. 

In order to check this scenario, we have performed nu- 
merical simulations of the same 2D binary mixture de- 
scribed above. We note that, on approaching a plastic 
event, caution must be taken to correctly minimize the 
energy without using a quadratic approximation. We 
observe that (i) the mode \&* is localized close to failure 
while (ii) 3(7) — ► S(7 C ) remains noisy and weakly corre- 
lated with the normal modes. As a consequence of this 
observation, £*('y) has a random, but typically non-zero 
limit when 7 — > j c . The non-affinc field is dominated by 
— £*(7c)/A*(7) x **(7 C ) and the non-afnne correction to 
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FIG. 3: Left: relative participation of the lowest nor- 
mal mode in the non-afnne elastic displacement field, a* = 
(£*/A*) 2 / 53 (£p/<V) 2 (dotted); lowest eigenvalue of the dy- 
namical matrix (solid); next several eigenvalues (dashed), b) 
In log-log scale (as a guide to the eye, the thick black line 
is V7= ~~ 7) : I/A 4 (circles); lowest eigenvalue (squares); next 
several eigenvalues minus their terminal values (diamonds). 



elasticity by /x ~ — (£,*) 2 /\* which diverges toward —00 
(see figure^]). In contrast, the Born term-which depends 
on the pair-correlation only- does not present any diver- 
gence. Since /x = /i Born + /x, on approaching failure, the 
system reaches a point 70 < 7 C at which Ji — — /x Born , 
whence /i vanishes. For 7 £ [7oj7c] 7 the shear stress is a 
decreasing function of 7: the material is unstable to any 
constant applied stress. This region is accessible to us 
because deformation-and not stress-is prescribed. 

In order to understand more specifically how the elas- 
tic constants behave close to j c , let us consider the 
functions Apf/y), which are continuous on a small in- 
terval close to 7 C (the second derivatives of the poten- 
tial are supposed to be regular). Close to the yield 
point, 7c, the deformation is dominated by the low- 
est normal mode: r(7) — r(7 c ) ~ #(7) (7J . (We 
project the deformation on the mode *ff* at the yield 
point.) From this relation and lp|l. we obtain the dom- 
inant contribution: dx/dj ~ — £*(7c)/A*(7). The coor- 
dinate x measures a true displacement in configuration 
space: we expect that no singular behavior occurs in this 
rescaled coordinate whence, \*(x) should vanish regu- 
larly, A* (a;) ~ ax close to x = 0. From this assump- 
tion, it results: 2^(7) 



V 2 Z*hc)(lc-j)/a. This rela- 
tion controls entirely the behavior of all observables when 
approaching the yield point: any observable A which be- 
haves regularly as a function of x (any regular function 
of molecular configurations) "accelerates" close to the 
yield point: dA/d'y ~ 1 /^j c — 7. In particular, we ob- 
tain, dr/d 7 ~ # P / % /7c Tr 7, and A* (7 ) = ^2a£*(~/ c - j), 
whence, ju ~ — (^*) 3 / 2 / ^J2a(^ c — 7). We could observe 
these scalings numerically by a careful approach to the 
yield point, as shown figure^ A similar divergence is ob- 
served for the compression modulus, but with a different 
prefactor, determined by the normal mode decomposition 
of the H field associated with pure compression. 

We now turn to the overall plastic event following fail- 
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FIG. 4: Left: Non-afRne elastic displacement field at a dis- 
tance 7c — 7 ~ 10~ 10 from the transition. Note the quadrupo- 
lar alignment with the direction of applied strain. Right: The 
local relative displacement field (the displacement of each par- 
ticle measured with respect to the average displacement of its 
neighbors) which is incurred after the entire plastic cascade. 
The solid line is a guide to the eye oriented along the oblique 
Bravais axis. 
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FIG. 5: Evolution of the displacement field during the irre- 
versible cascade corresponding to the event circled in figure^ 
Top: contribution a* (solid line) of the critical mode and a p (t) 
(dashed line) of the next five modes to the displacement field. 
Bottom: sum of the squares of the forces on the particles. 

ure (see figure 0] and [SJ) . We have already shown, in 
similar atomistic systems, that any single plastic event 
involves a cascade of local rearrangements. |7j Our pre- 
ceding work suggested that the overall plastic event was 
controlled by long range elastic interactions and differed 
from the displacement fields which dominate the onset 
of failure. Our present normal mode decomposition al- 
lows us to gain more insight into this process. Writing 
r(t) — f (0) = J2 p A£ p (t) * p (7c), we extract the quantities 
a p = (A£ p (t)) 2 / £ p (A£ p (i)) 2 . which are shown figureEl 
for the lowest frequency modes. To trigger the relax- 
ation, we shear the system forward by a small amount 
of shear, 7 — j c ~ 10~ 5 . The initial affine displacement 
serves as a perturbation and projects randomly on the 
normal modes, whence the contributions a p start around 



zero. We observe that (i) the initiation of the cascade is 
clearly dominated by the critical mode (ii) this effect sud- 
denly stops before reaching the first peak in £\ Ff (this 
peak correspond to the first inflection point of energy vs. 
minimization step) (iii) the subsequent displacement ap- 
pears to be random, when projected on the lowest part 
of the spectrum, indicating that low frequency normal 
modes are irrelevant for the latter stages of plastic fail- 
ure. The system escapes its initial inherent structure in 
the direction of the lowest normal mode, whereas, in the 
late stages of a plastic event the emergence of long-range 
elastic interactions supersedes the low energy modes re- 
sponsible for the onset of failure. 

To conclude, we wish to stress that we expect our anal- 
ysis to apply to many materials, to modes of deforma- 
tion other than uniform shear, and to several experimen- 
tal protocols at the nanoscale, including nanoindentation 
studies. Knowing the detailed behavior of elastic con- 
stants at incipient plasticity opens the route toward a 
possible control of material deformation at the nanome- 
ter scale. 
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